(N 
O 



Abstract. There exists both theoretical and observa- 
tional evidence that the magnetic field decay in neutron 
stars may proceed in a pronounced non-linear way during 
a certain episode of the neutron star's life. In the pres- 
ence of a strong magnetic field the Hall-drift dominates 
the field evolution in the crust and/or the superfluid core 
of neutron stars. Analysing observations of P and P for 
sufficiently old isolated pulsars we gain strong hints for a 
significantly non-linear magnetic field decay. Under cer- 
tain conditions with respect to the geometry and strength 
of a large-scale magnetic background field an instability 
is shown to occur which rapidly raises small-scale mag- 
| netic field modes. Their growth rates increase with the 
£N) ■ background field strength and may reach ~ 10 4 times 
the ohmic decay rate. Consequences for the rotational and 
^ ■ thermal evolution as well as for the cracking of the crust 
of neutron stars are discussed. 
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1. Introduction 

There is still an ongoing scientific debate whether the mag- 
netic fields of pulsars decay at all (see, e.g., Elcgimbau & 

and references therein). Thus 



de Freitas Pacheco 2001 



the question whether this questionable decay proceeds in 
a non-linear way seems to be somewhat academic. How- 
ever, the neutron star magnetic fields are the strongest 
known in the universe and already this very fact suggests 
the idea that, if there is a field evolution at all, it should 
be a non-linear one. As we will show, there exist both the- 
oretical and observational evidences that the field decay 
proceeds essentially in a non-linear way, at least in certain 
regions of the star and surely only during certain episodes 
of its life. 

In order to be able to deal with well defined phys- 
ical conditions we confine ourselves here to considering 



— therefore, there is a stable compactness ratio M/R and 
no screening of the magnetic field can happen. 

All these constraints are met by the majority of neutron 
stars the age of which is > 10 5 years. Therefore, if any 
thermal and/or rotational evolution of SOINSs is observed 
beyond that predicted when assuming a constant magnetic 
field, it is with a high probability caused by an evolving 
one. 

Here, we intend to consider as the non-linearity of the 
field evolution the Hall-effect or Hall-drift, as it occurs in 



the crust (see, e.g., Shalybkov & Urpin 1997), i.e., both 



ambipolar diffusion and any convcctive motion of the neu- 
tron star matter is excluded. 

The effect of the Hall-drift on the magnetic field evolu- 
tion of isolated neutron stars has been considered by a 
number of authors (see, e.g., Haensel et al. 1990; Gol 



jlution of sufficiently old isolated neutron stars drcich fc Reisencggcr 1992]; [Muslimov 1994|; [Naito fc Ko 



the ev 

(SOIN is). We characterize them by the following prop- jima 1994|; |Shalybkov fc Urpin 1995|, |l997j ; |Urpin fc Shaly 
erties: | bkov 19951 , EUj i |Vainshtcin et al. 2000|) . They discussed 



the supernova fall-back accretion phase is finished and 
a stable density stratification p(r) has been estab- 
lished; 

the rediffusion of the magnetic field has been com- 
pleted; 



the redistribution of magnetic energy from an initially 
large-scale (e.g., dipolar) field into small-scale compo- 
nents via the Hall-term. Though the Hall-drift itself is a 
non-dissipative process, these changes in the field geom- 
etry may in principle accelerate the field decay However, 



- the rotation is slow enough so that the irradiation of (1994), ghafybkov fc Urpin (1997), and Urpin fc Shaly 



the re su lts presented in Naito fc Kojima (1994 ^ Muslimov 



gra vitational waves docs not contribute to the spin 
down; 

the appearance of glitches is less probable; 



bkov (1999 ) suggest the conclusion that the field decay is 
not modified drastically. 

Goldrcich & Reisencggcr Goldrcich & Rcisenegger 



the crust is almost completely crystallized; 



ide, 



when 



the temperature profile has become flat enough, to 
avert convection in the outermost thin liquid layer and 
in the core region as well as the occurrence of a ther- 
moelectric instability in the crust; 
the electron relaxation time is already so large that 
the magnetization parameter may reach values large 
enough to make the electron transport processes non- 
linear: 



(1992) developed the idea of the Hall-i 
starting with a large-scale magnetic field small-scale field 
components are generated down to a scalelength ^ cr it, 
where the ohmic dissipation begins to dominate the Hall- 
drift. 

In some of the above-mentioned investigations numer- 
ical instabilities arc mentioned if cither the field structure 



— the absence of accretion, as a consequence of isolation, 



ensures that there are neither external sources of an- 
gular momentum nor of heating; 



(Urpin & Shalybkov 1995|) or the 


g( 


Naito & Kojima 1994 


; Urpin & 



Shalybkov 1999). Also, when considering the thcrmomag- 
nctic field generation in th e crust of young neutron stars 
( Wiebicke fc Geppert 1996 ), where small-scale modes are 
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the first ones to be excited, numerical instabilities oc- 
curred exclusively caused by the Hall-drift. 



Recently we have shown (Rheinhardt & Geppert 2002) 
that all the observed instabilities are in their essence very 
likely not of numerical origin but have physical reasons: 
A sufficiently strong and inhomogeneous large-scale back- 
ground field is unstable with respect to small-scale per- 
turbations, which rise very rapidly in comparison with 
the relevant ohmic decay time of the background field. 
This rapid transfer of energy from the background into 
small-scale field components proceeds not via a cascade 
but jump-like across wide spectral distances. The unsta- 
ble perturbations show small radial structures close to the 
neutron star's surface whereas their lateral structures are 
of medium size. Of course, this Hall-instability can act 
only during an episode of the field decay, which unavoid- 
ably leads to a zero field. This episode, however, may have 
observable consequences. 

In the following we will present the theoretical and ob- 
servational hints for the non-linear field decay in isolated 
neutron stars. After a description of the Hall-drift induced 
field instability we discuss its effects on observable quan- 
tities of the above introduced SOINSs. 

2. Theoretical motivation 



Under which conditions becomes the evolution of the mag- 
netic field non-linear? The field evolution is determined by 
Maxwell's equations in quasi-stationary approximation 



B = -c curl E , divB = 



47T 



curlJS (1) 



together with Ohm's law, which has in the absence of con- 
vectivc motions the form 



E 



Rj 



(2) 



In the presence of a magnetic field the electric conduc- 
tivity, <7, and the resistivity, R, respectively, become ten- 
sors whose components parallel and perpendicular to the 
magnetic field along with their Hall-component are rep- 
resented by 




R- 



R± Rh 
- i?H R±. 
R\\ 



> (3) 



if B goes along the z-axis. 

Let us consider the evolution of a magnetic field com- 
pletely confined in the crust of the neutron star. There, 
the field could have been generated either by convection 



in the outermost layers of the proto-neutron star or by 



a thermoelectric instability during the early period of the 
neutron star's life, when the temperature gradient in the 
liquid crust is extremely large (see, e.g., Urpin ct al. 1986; 



Wiebickc & Geppert 1996). A peculiarity of the transport 



processes in the crust is that the components of the resis- 
tivity tensor parallel and perpendicular to the magnetic 



field coincide (see Yakovlev & Shalybkov 1991), because 
the electrons move as a degenerate gas through the crystal 
lattice of the crust. In relaxation time approximation they 
found 



(Tq : 



00 



1 + (ujbt) 2 



ch = lubt o j 



(4) 



where lob — eB /m* e c is the Larmor frequency of the elec- 
trons, m* their effective mass, and r their relaxation time 
between two collisions with the phonons or impurities of 
the crustal lattice. As a consequence of Eq. (|J) the paral- 
lel and perpendicular components of the resistivity tensor 
just coincide 



Rn = Ra 



Rn 



oh 



(5) 



The scalar electric conductivity o~o as well as the relaxation 
time r are complicated functions of the density p(r), the 
temperature T(r, i), the impurity concentration Q and of 
the chemical composition (i. e., charge number Z( r) and 
mass number Air) 



see, e.g. 



[Jrpin fc Yakovlev 1980) 



Ohm's law (|2|) together with Eq. (|3|) allows of repre- 
senting the electric field by use of the unit vector of the 
magnetic field es = B/\B\ as 

E = R\\{e B ■ j)e B +i?±(e B x j) x e B + Rn(j x e B ),(6) 
which because of R± = R\\ yields 

E = Rp + R H (j x e B ) = (Jn 1 (j + uJ B r(j xe B )). (7) 

Taking now the curl of E and using Ampere's law we 
arrive at the Hall-induction equation for the crystallized 
neutron star crust: 



B = curl 

47T 



curlJB +uj B T{cm\B Xe B )) (8) 



which makes immediately clear that the field decay will 
proceed in a significantly non-linear manner as soon as 
(and where) 



eB(r,t) 
m*(r)c 



(p(r),T(r,t),Q,A(r),Z(r)) > 1 , (9) 



i.e., the significance of the non-linearity depends not only 
on the field but on all the complex physical conditions in 
the crust. 

In order to get an impression how the magnetization 
parameter ojbt in the crust of an isolated neutron star 



evolves we present here results obtained by Page et al 
(2000D for a neutron star whose state of matter in the 



core is described by a medium equation of state. The ini- 
tial polar surface field strength is assumed to be 10 13 G, 
the density up to which the field is initially penetrating is 
po = 10 14 g cm" 3 ; for the impurity concentration Q = 0.01 
is supposed and the chemical composition of the crust 
should be that of cold catalyzed matter. Since the field de- 
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Fig. 1. The magnetization parameter wgr as a function 
of the density at different moments of neutron star's life. 



cay as presented by Page et al. (2000|) has been calculated 
just by neglecting the non-linear effects, i.e., considering 
only the linear ohmic decay, Fig. |l| does not present the 
"real" evolution of ujbt because B(r,t) certainly evolves 
differently under the influence of a strong Hall-drift. Nev- 
ertheless, these results may serve as a strong hint that the 
latter may play an important role for the magnetic field 
decay in neutron stars. When, e.g., lubt becomes as large 
as 1000 as a consequence of the star's cooling down due 
to which the electron relaxation time rises by orders of 
magnitude, a linear description of the field decay is surely 
no longer justified. 

3. Observational motivation 

The best known observational quantities of SOINSs are 
obtained when they appear as radiopulsars: their pulsa- 
tion period P = 27r/0 (O being the related rotational 
velocity), its time derivative P and, unfortunately only in 
a few cases, its second time derivative P. The spin-down 
behaviour is to a high accuracy determined by the require- 
ment, that the loss of rotational energy E lot = I0 2 /2 (I is 
the moment of inertia of the neutron star) is equal to the 
power of the magneto-multipole radiation E mmr , emitted 
by the rotating magnetized neutron star appearing as a 
radio pulsar: 



E IC 



-EL 



(10) 



As already discussed above, for pulsars older than w 10 5 
years, the spin-down as a consequence of the emission 
of gravitational waves can surely be neglected, too. The 
loss of angular momentum due to a relativistic wind of 
charged particles (caused by the pulsar mechanism) shows 
the same functional dependences on the neutron star ra- 
dius R, the magnetic field and the rotational velocity as 



the magneto-multipole radiation (Shapiro & Tcukolsky 
1983). Clearly, Eq. (|l0|) is not valid during glitches. 



According to Krolik (1991 ) the power of the magneto- 
multipole radiation is given by 

2(1+1) 



E u 



= —ys 

l,m 



2 
I in 



2H\ 

w 



R 2 (\B lm \y 



(11) 



where (\Bi m \) is the average of the multipole component 
of the magnetic field, defined by its multipolarity I and 
projection m, across the spherical surface. Si m is a con- 
stant reflecting the geometry of Bi m (sec Krolik, Krolik 
199ip. S ince R and I are fixed for SOINSs, we have with 
Eqs. |10| and ( pd] ) an unambiguous relation between the ob- 
servational quantities P and P and the magnetic field. In 
other words, if we were able to calculate all the multipole 
components of the magnetic field as functions of time and 
if these together with the measured quantities P(t) and 
P(t) satisfied Eq. (10), we had understood the field de- 
cay. However, as a consequence of the sum occurring in 
Eq. ([ll]) it is impossible to infer from P(t) and P(t) all 
the (|I?i m |) in an unambiguous way. Therefore, we con- 
sider here only one term of the sum to be responsible for 
the spin-down and denote the corresponding (|Bj m |) sim- 
ply by B. (Usually, the dipolar term, I = m = 1, is chosen.) 
Then, by equating Eq. ( |To| ) with Eq. ([TT|) , the change of 
the neutron star's rotation can be written as 



O = -KQ r 



(12) 



where the constant K comprises the geometrical constant 
Sim as well as radius and moment of inertia and the mag- 
netic field. The exponent n = 21 + 1 is called the "braking 
index" , where n = 3 for a pure rotating dipole. The stan- 
dard way to calculate the braking index is to take the time 
derivative of Eq. (O) resulting in 



_ OO _ PP 
" ~ 02" ~ ~ P2" 



(13) 



However, O (or P) can be obtained with sufficient accu- 
racy for a few young pulsars only. 

Johnston & Galloway (1999) suggested an alternative 
way to calculate n without employing P. By direct inte- 
gration of Eq. (|l2|) over a time interval T 3> P 
,n 2 m 

-KT (14) 



dO 

O" 



they calculate the braking index by 



71 = 1 



Oi0 2 T 



(15) 
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i.e., knowing the rotational period and its time derivative 
of two observations, separated by T, one is able to ob- 
tain the braking index without needing P. For the Crab 
pulsar and B 1509-58 they found a good coincidence of 
the values of n obtained with both methods. Thus, there 
exists a certain warranty that the application to other ra- 
dio pulsars will yield reliable results, too. In their Tabic 1 



Johnston & Galloway (1999) present the braking index n, 
calculated according to Eq. (|l5|) along with P and P for 20 
pulsars. Six of them have negative braking indices which 
possibly could be understood as due to the occurrence of 
glitches or the continuation of a post-glitch recovery dur- 
ing the observational periods. As will become clear later, 
a decay of the NS magnetic field can never be responsible 
for these negative values (whereas a field growth can) . The 
eight pulsars with the smallest error bars for n (B0540+23, 
B0611+22, B0656+14, B0740-28, B1915+13, B2002+31, 
B2148+52 and B2334+61) have all positive braking in- 
dices, five of them with n > 20. Because of their simi- 
lar braking indices B0919+06, B1221-63 and B1907+10 
may belong to that group too, but their error bars are 
much larger. Two further pulsars, B0154+61 and B1356- 
60, with n > 3 have unacceptably large error bars. By use 
of a recently available pulsar catalogue (Manchester et al. 
200l|y]ve found two more SOINSs (B0727- 18, B1822-09) 
for 

1993|)7]'l'heir derived values ot n are the largest tound by 
the above method, although the error of n for B1822-09 is 
rather high (about 60%). All pulsars mentioned have an 
active age r a = P/2P between 4.13 ■ 10 4 (B2334+61) and 
1.7 • 10 6 (B1907+10) years, i.e., are just at the range of 
ages, the SOINSs should be. For a collection of the selected 
SOINSs' data see Table |. 

A basic supposition for both methods to calculate the 
braking index is that the coefficient K in Eq. (|l2|), by 
virtue of Eq. (11) proportional to B 2 , is really a constant. 



justification for that assumption can only be obtained for 
each individual pulsar by considering its actual time series 
of P and P. Re-evaluating now Eq. (|IJ|) with n replaced 

by n we find with P = JCBP- n+1 (2BP - (n - 2)BP) 
the correction of n to the "observed" n to be proportional 
to the time derivative of the magnetic field: 

n = n - 2^- = n - 2 sign(SP) — . (18) 
BP t b 

Here wc introduced rp and tb as characteristic times of 
the rotational period and the magnetic field, respectively, 
but we discuss this relation for a spin-down situation only, 
that is, P > How should this equation be interpreted? 
In case there is no magnetic field change, that is B = 0, 
pulsars with n < 3 are not understandable in the frame- 
work of the magneto-multipolc radiation model, because 
the neutron star will certainly not be a rotating monopole 
(with small higher multipole contributions). In case we as- 
sume that the spin-down is completely due to magneto- 
dipole radiation, an assumption perhaps not completely 
wrong for SOINSs, we have n = 3 and n < 3 indicates an 
increasing field while n > 3 is the signature of a decaying 
one. The interpretation of n < 3 as caused by an increas- 
ing field is supported by the fact, that such braking indices 



are observed especially for very young pulsars ( Lync et al 



ch also earlier observations exist (|laylor et al.| 1996|), in which both rediffusion of a field submerged din- 



ing the supernova fallback (Geppert et al. 1999) and/or 
thermoelectric field generation may still proceed. Alterna- 
tively, deviations of n from its classical value 3, even values 
of either sign up to 10 6 , can be understood as due to in- 
ternal frictional instabilities occurring between the crust 
and the supcrfluid, almost independently on the evolu- 



tion of the NS magnetic field ( Shibazaki fc Hirano 1995 ; 
IShibazaki fc Mochizuki 1995| ; |Kimdt 2001|) . However, this 



Now let us relax that assumption in so far as we will allow 
for a varying field while the other quantities entering K re- 
main constant, in agreement with our concept of SOINSs. 
Then we can rewrite Eq. ( |l2|) into 



P = KB 2 {t)P 



-n+2 



(16) 



where for a standard neutron star (R = 10 6 cm, / = 
10 45 gcm 2 ) the constant K « 10~ 39 cm s 3 g" 1 if n = 3. 
We will nevertheless assume in the following that the ro- 
tational evolution is during a period T still quite well 
described by a power-law similar to Eq. ([12]), where, of 
course ,. the quantities K and n have to bo replaced by 
modifi ed ones: 



explanation applies, except for strong vortex pinning, only 
for older (age > 1.8 10 7 yrs) NSs. Thus the assumption of 
a rapid field decay could just explain large observed values 
of n for younger objects, as the SOINSs listed in Table [l] 
are. 

Also processes in the magnetosphere may cause values 
of fi different from 3, but these effects are considered to 
cause h < 3 in young pulsars (see, e.g., Melatos 1997; 
|Casini fc Montcmayor 1998 ) 



In this context it should be mentioned that the con- 



clusion of Rcgimbau & dc Frcitas Pachcco (2001) a decay- 
ing magnetic field causing a braking index greater than 
three was "in disagreement with observations" is obvi- 



ously wrong. Along with those for which Johnston & Gal 



loway (1999| ) calculated the braking index it has been dc 



(17) 



The description of the rotational behaviour by a power law 
should be sufficiently precise as long as the time T (see 
Eq. ([l5|)) is small in comparison with the characteristic 
time of the magnetic field evolution. Of course, a final 



termined for only five very young pulsars employing P. In 
all there exists only a single case (J1119-6127) with almost 
exactly n = 3. Thus, the assumption of a varying (either 
increasing or decreasing) magnetic field is in general the 
more appropriate one. 

Here we intend to investigate, whether the observed 
relatively large values of n may indicate a non-linear field 
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Table 1. SOINSs supposed to be in the phase of non- 
linear magnetic field decay. Rotational period P and its 



time derivative P are mainly taken from Table 1 of John- 
ston fc Galloway (1999| ) 



n is the rounded mean value 
of n given there. The data for B0727-18 and B1822-09 



come from the pulsar catalogues Taylor et al. (1993) and 
Manchester ct al. (2001). The active age, polar dipole field 
strength and field decay rate are calculated as described 
in Sect. 



PSR(B) 


P 


P 


h 




B d 


B d 




M 


[io- 14 ] 




[10 s yrs] 


[lO^G] [10 s G/yr] 


0540+23 


0.246 


1.531 


12 


2.55 


3.93 


0.35 


0611+22 


0.334 


5.921 


20 


0.89 


9.00 


4.28 


0656+14 


0.385 


5.507 


15 


1.11 


9.32 


2.52 


0727-18 


0.510 


19.01 


319 


4.26 


6.30 


11.66 


0740-28 


0.167 


1.687 


26 


1.57 


3.40 


1.25 


0919+06 


0.431 


1.375 


29 


4.97 


4.93 


0.65 


1221-63 


0.216 


0.493 


19 


6.94 


2.09 


0.12 


1822-09 


0.769 


52.43 


142 


2.33 


12.8 


19.20 


1907+10 


0.284 


0.265 


24 


17.6 


1.75 


0.05 


1915+13 


0.195 


0.724 


36 


4.28 


2.40 


0.46 


2002+31 


2.128 


7.562 


23 


4.46 


25.7 


2.88 


2148+52 


0.332 


1.005 


50 


5.24 


3.70 


0.83 


2334+61 


0.495 


19.02 


9 


0.41 


19.6 


7.13 



decay. For that purpose let us consider as a typical ex- 



ample B0740-28, found in Table 1 of [Johnston & 
loway 



1999) 



Gal- 



. with n w 26. From its observed P = 0.167s, 
P = 1.68 • 10~ 14 we can infer a dipolar surface magnetic 
field strength B d = 6.4 • 10 19 y/p~P « 3.4 • 10 12 G which 
is in the range of values typical for such "middle aged" 
radio pulsars. Using the standard formula for the active 
age T a = P/2P we find T a w 1.6 ■ 10 5 years, but this relies 
on the assumption that the rotational evolution follows 
all the time the power law ( |i~7| ) with n = 3. Of course, 
this value of r a is at best an upper limit for the star's age, 
since we just argue in favour of an n > 3 during some (un- 
known) period of the rotational evolution still lasting. On 
the other hand a much shorter active age is unlike because 
none but one of the mentioned pulsars is reported to be 
associated with a supernova remnant. The only exception, 
B2334+61, with r a = 4.1 • 10 4 yrs, is associated with the 
fairly old SNR G114.3+0.3 |Bccker et al. (1996|); thus the 



active age seems to be a sound estimate of the real age. 

If we exclude the possibility of magnetic field decay 
we have for our example (B0740-28) n = n«26, i.e., the 
neutron star's magnetic field structure must be that of a 
2 12 -pole because of I fs 12. As a consequence of Eq. ( p+| ) 
this implies, that the surface magnetic field strength of 
such a field must be many orders of magnitudes stronger 
than that of a dipolar field in order to explain the ob- 
served spin-down. Therefore, the assumption of a constant 
magnetic field seems not to be very likely. If we instead 
assume that the spin-down is caused by magneto-dipole 
radiation, i.e. n = 3, we can easily estimate the rate of the 



magnetic field decay. Rearranging Eq. (|l^) and inserting 
the observed or inferred values for P, P and we find 



Ba 



B d P 
2P 



4G/s« 1.25- 10 8 G/yr. (19) 



This rate of the dipolar-field decay should now be com- 
pared with the typical ohmic decay rate of a crustal mag- 
netic field. The ohmic decay time of a crust of thick- 
ness L = 10 cm and an averaged electric conductivity 
cr = 5 • lO^s" 1 is 



To hii 



2.2 • 10 9 yrs, 



(20) 



providing an upper bound for the real field decay times. 
By use of the above estimated B d = 3.4 ■ 10 12 G we find 



B d = 1.25 • 10 8 G/yr > 



1.5-10 3 G/yr. (21) 



This comparison is robust, i.e., even if we somewhat over- 
estimated L or <7o (which, because of the unknown impu- 
rity concentration, could indeed also be underestimated), 
the magnetic field decay proceeds obviously much faster 
than would be predicted for a (linear) ohmic decay (see, 
e.g., Tauris & Konar 2001). What has been demonstrated 



here for a single example proves to be valid for the com- 
plete selection of SOINSs specified above, see Table § 

As demonstrated in Sect. |^ the only possible modifica- 
tion of the ohmic decay in the context of our model is due 
to the Hall-drift. That is, any significant acceleration of 
the magnetic field decay is explainable then and only then, 
if the Hall-drift can be shown to cause time scales of the 
field evolution significantly smaller then the ohmic one. In 
order to visualize the clear signs of non-linear field de- 
cay for the selected SOINSs we present in the upper panel 
of Fig. H the correlation of the characteristic time of the 
magnetic field tb = \B d /B d \ with the magnetic field B& 
itself. Bearing in mind, that thirteen objects form only a 
very small statistical ensemble, there nevertheless seems 
to exist a pronounced trend indicating non-linearity of 
the field decay: Evidently, tb falls significantly with in- 
creasing B d which is compatible with the idea that the 
stronger the field the faster proceeds the redistribution of 
energy into smaller modes by the Hall-drift and, conse- 
quently, the faster the energy is drained from the dipolar 
field mode. Note, that even the two pulsars (B0727-18, 
B1822-09) with apparently extraordinary high values of 
B d fit well into this trend. 

In addition, we show in the lower panel of Fig. ^| the 
correlation of tb with the active age t q suggesting the 
conclusion that the acceleration of the field decay by the 
Hall-drift acts only temporarily. This makes sense because 
that effect of course destroys its own basis, namely the 
dominance of the Hall-drift over the ohmic decay. Cer- 
tainly it is to premature to take the link between the field 
decay and the observations as a proven fact. However, the 
apparent trend gives enough reason to pursue this idea. 
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4.0X1CT r 



noted, this analogy is incomplete at least as far as one curl 
operator has to be replaced by its inverse, which makes the 
conclusions about the Kolmogoroff spectrum questionable. 

Attempts performed by the present authors to solve 
the non-linear Eq. (||) in a full sphere assuming a con- 
stant density profile reproduced the results of Shalybkov 



& Urpin (1997) for moderate values of the magnetiza- 
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Fig. 2. Characteristic time of the magnetic field, tb, ver- 
sus the magnetic field B itself (upper panel) and versus 
the active age r a (lower panel) for the SOINSs given in 
Table [j] 

In the following chapter we will present an instabil- 
ity which proceeds on time-scales much smaller than the 
ohmic scale. 



4. A Hall drift induced magnetic field instability 

4-1- Hints for the occurrence of an instability 

Now, let us return to Eq. (|8|) which together with div B = 
and appropriate boundary conditions determines the 
magnetic field evolution in the crust of isolated neutron 
stars. It is generally accepted that the non-linear Hall- 
term causes a re-distribution of magnetic energy, when 
it is at a certain "initial" state stored in the mode char- 
acterized by the largest spatial scale, say a dipolar field 



mode. Goldreich & Rciscncggcr (1992) coined the term 



Hall-cascade for that process by stressing an analogy to 
the vorticity equation for an incompressible fluid. As they 



tion parameter u>bt. Starting with a dipolar field modes 
of smaller scales were excited and hclicoidal waves were 
observed. But, as soon as ujbt exceeded a certain value, 
typically about 100, the spectral resolution, attainable by 
the numerical code was no longer sufficient to ensure con- 
vergence because the amplitudes of small-scale modes at 
the tail of the resolvable spectrum started to grow very 
quickly. This motivated us to search for possible physical 
reasons for that apparent numerical instability. 

For simplicity we assume the conductive properties of 
the matter to be constant in space and time, that is, 
we assume constant Co and r e /m*. Thus, Eq. (|^) can 
be rewritten in dimcnsionless variables such that it no 
longer contains any parameter and the evolution of the 
magnetic field is solely determined by its initial configu- 
ration B(x, 0). For that purpose we normalize the spatial 
coordinates by a characteristic length L of the model (e.g., 
the neutron star's radius or the thickness of its crust), the 
time by the ohmic decay time 4iraoL 2 /c 2 and the mag- 
netic field by = m* e c/eT e . The governing equations in 
these dimensionless variables read 

B = AB - curl(curlB x B ) , divB = , (22) 

where the differential operations have to be performed 
with respect to the now dimcnsionless spatial and time 
coordinates x, y, z and r, respectively. 

Let us first have a look onto energetics. Using standard 
arguments one can immediately state that in the absence 
of currents at infinity the total energy of any solution of 
Eq. (^U) is bound to decrease monotonically to zero since 
the Hall-term curl B x B is unable to deliver energy (nor 
to consume it). That is, all these solutions are asymptot- 
ically stable in a strict mathematical sense following Lya- 
punov's definition of stability. However, this applies also 
to a lot of slowly time varying processes (background) for 
which nevertheless the transient occurrence of instabilities 
is a generally accepted idea (e.g., the temporary convec- 
tion in a cooling liquid). Therefore, it seems reasonable 
to adopt as a criterion of instability the following: Define 
a steady reference (equilibrium) state in the spirit of the 
usual stability analysis to be a snapshot of the background 
process. If the analysis detects growing, that is, unsta- 
ble perturbations with a time increment large enough in 
comparison with the time decrement of the background 
process under consideration the latter may be accounted 
unstable in a short enough period of time around the in- 
stant of the snapshot. 

Of course, when adopting this concept a necessary con- 
dition for the occurrence of an instability is as well the 



U. Geppert & M. Rheinhardt: Non-linear Magnetic Field Decay 



existence of terms in the linearized equation which are 
potentially able to deliver energy to the perturbations. 
Linearization of Eq. (22) about the background field Bq 
yields 



z = +l 



b = Ab — curl( curl B x b + curl b x B ) 
div b = 



(23) 



describing the behavior of small perturbations b of the 
reference state Bq. Indeed, here along with the energy- 
conserving term curl b x Bq now as a second Hall-term 
curl Bq x b occurs which may well deliver or consume en- 
ergy (from/to b) since in general the integral J v (cm\ Bq x 
b) ■ cmlbdV will not vanish. Of course, this is due to the 
fact that the linearized equation describes the behavior of 
only a part of the total magnetic field. Actually, pertur- 
bations may grow only on expense of the energy stored in 
the background field. 

Considering Eq. (|23|), we can determine a scale be- 
low which the ohmic dissipation dominates the Hall-drift. 
Estimating | curl Bo | and |curlb| by B and b/l, respec- 
tively, we find the critical scale of b to be Z cr ;t < 1/Bo, 
which is de-normalized / cr jt < L/(ujg n T e ), identical with 
the expression derived by Goldreich & Rcisencgger (1992) 



considering the Hall-cascade in analogy with the turbu- 
lent flow of an incompressible fluid. 

4-2. Analogy with the Vt — (u> x j)~dynamo 

A strong support for the idea of a Hall-drift induced mag- 
netic field instability is provided by the analogy of the 
linearized Hall-induction equation (^3|) to the mean-field 
induction equation including the so-called (u> x j)-cffcct 
found by [Radler (1969)): 



B = AB curl ( -u n x B + curl Bxu). (24) 

Here, B is the averaged magnetic field, uq describes a 
(large-scale) differential rotation and a; is a turbulence- 
related coefficient, which reflects the offset of the isotropy 
of the turbulence by Coriolis forces. Eq. ( jg3 ) has been 
proved to allow for magnetic instabilities (i.e., dynamos; 



see |Rajilcr 1969|T 



Comparing Eq. (23) with Eq. (|24|) when letting corre- 
spond b to B it becomes clear, that curl Bo (then corre- 
sponding to —uq) must not be a constant, since —uq is 
a shear velocity. This implies that not all second deriva- 
tives of Bo with respect to the spatial co-ordinates may 
vanish. It should be noted that the analogy lacks com- 
pleteness only insofar that the dynamo acts even with a 
constant u>. Note that u> enters the term which is energy 
conserving. 

4-3. Mathematical treatment 

Let us now specify the geometry of our model and the 
background field. We consider a slab which is infinitely 




z = -1 

Fig. 3. Sketch of the model and background field geom- 
etry. In z < — 1 a perfect conductor, in z > 1 vacuum is 
assumed. 

extended both into the x- and y directions but has a finite 
thickness 2L in z-dircction. 

The background field Bq is assumed to be everywhere 
parallel to the surface of the slab pointing, say, in x— 
direction and to depend on the z-coordinate only, i.e., 
Bo = f{ z ) e x {s-x being the unit vector in x— direction) . 
Note, that by this choice curl B x Bq represents a gradi- 
ent. Thus, the unperturbed evolution of the background 
field is not at all affected by the Hall-drift; in the absence 
of an e.m.f. it would decay purely ohmically. 

Further on we decompose a perturbation b into its 
poloidal and a toroidal components, b — bp + bx, which 
can be represented by scalar functions S and T, respec- 
tively, by virtue of the definitions 

bp = - curl ( e z x VS) , b T = -e z x VT , (25) 

ensuring div b = for arbitrary S, T. 

For the sake of simplicity we will confine ourselves to 
the study of plane wave solutions with respect to the x- 
and y-directions, thus making the ansatz 



(x,t) 



t 



(z) exp (ikx + pr) , (26) 



where k = (k x ,k y ) ^ 0, x = (x,y) and p is a com- 
plex time increment. This ansatz guarantees as well the 
uniqueness of the poloidal-toroidal decomposition since 
from A (S, T) = it follows (S,T) =0 with A being the 



two-dimensional lateral Laplacian (see Briiuer & Radlei 
1988| ). With Eq. @ we obtain from Eq. (pj) two coupled 



ordinary differential equations 

pt - t" + k 2 t = \k x f{s" - k 2 s) - ik x f"s ^ 

ps — s" + k 2 s — ikyf's = —ik x ft , 

where dashes denote derivatives with respect to z. To- 
gether with appropriate boundary conditions Eqs. (27) 
define an eigenvalue problem with respect to p. 

Having in mind that the neutron star's crust is neigh- 
boured upon the superconducting core on the one and a 
region with very low conductivity on the other side we 
choose the following simplified boundary conditions: Out- 
side the star we assume vacuum, that is, curl B = above 
the slab. Because of the finite conductivity inside the slab 
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Fig. 4. Growth rate as a function of k x and k y for Bq = 
1000. Negative values were set to zero. 

we require continuity of all components of B across this 
boundary. Considering the crust-core boundary an elec- 
tric field must be prevented from penetrating into the per- 
fectly conducting region below the slab, z < — 1, that is, 
the normal magnetic and tangential electric field compo- 
nents must vanish. In terms of the scalars s and t this 
means [s] = [s'] = t = for the vacuum boundary and 
s = if = for the perfect conductor boundary where [.] 
denotes the jump across a boundary. Note, that for t' = 
to be valid the vanishing of B at the perfect conduc- 
tor boundary is required. Making use of the vacuum solu- 
tion for the upper halfspace, z > 1, vanishing at infinity, 
the vacuum boundary condition for s can be expressed as 
s' = — ks at z = 1, with k = |fc|. 

Finally we specify the background field profile: Guided 
by the conditions under which the above discussed f2 — 
(u> x j')-dynamo may work we conclude that f(z) has to 
be at least quadratic in z. Taking into account the bound- 
ary conditions we choose f(z) = Bq(\ + z)(l — z 2 ). From 
now, the coefficient Bq is the only free parameter of the 
eigenvalue problem. 

4- 4- Numerical results 

For certain ranges of the wave numbers k x , k y and values 
of Bq > 3 we found eigenvalues p with a positive real part, 
i.e., exponentially growing perturbations. The dependence 
of the growth rate 5ft(p) on the wave numbers for Bq = 
1000 is shown in Fig. |[ 

Figure |5| shows the dependence of growth rate and 
wave number fc" lax of the fastest growing mode on Bq. 
An interesting feature is, that the maximum growth rates 
occur for all Bq considered at k y = 0. Of course this asym- 
metry is due to the choice of the background field: once 
it was chosen parallel to the y-direction the maximum 




Bo 

Fig. 5. Growth rate and wave number fc™ ax of the fastest 
growing mode as functions of Bq. Thick and thin lines 
correspond to growth rates and fc™ ax , respectively. 

growth rates would occur at k x = 0. Note that the most 
unstable eigenmodes are always non-oscillatory, though 
oscillating unstable ones exist. 

Evidently, for each Bq above a certain threshold (say, 
Bq = 50) there exists a range of wave numbers within 
which the obtained growth rates are in agreement with 
the constraint, formulated above: In comparison with the 
background field decay the growth of the most unstable 
perturbations is a fast process; thus we may consider it as 
"episodically unstable" . 

With respect to the asymptotic behavior uq — > oo 
for a fixed (unnormalized) background field one has to 
note that the time increment p is normalized on the 
ohmic decay rate (oc (Tq 1 ). From Fig. [5] it can be inferred 
3?(p) oc B^ , q < 1 for B > 100, which means that in the 
limit of negligible dissipation the growth rate in physical 
units tends to zero. 

Figure |^ shows the cigensolutions (s,t)(z) of the 
fastest growing mode for three different values of Bq. One 
can observe that with increasing Bq the toroidal field be- 
comes more and more small-scale and concentrated to- 
wards the vacuum boundary. In contrast, the correspond- 
ing poloidal field remains large-scale. The magnetic field 
structure of the fastest growing mode for Bq = 2000 is 
shown in Fig. |^. Assuming a ratio of crust thickness and 
neutron star radius R equal to 0.1 the shown lateral pe- 
riod of that mode fits ~ 20 times into the circumference 
of the star. 

5. Discussion and Conclusions 

Clearly, any assignment of the results gained by help of a 
very simplified model to real neutron stars has to be done 
with great care. Even when accepting the plane layer as a 
reasonable approximation of the crust one has to concede 
that the very specific profile of Bq assumed above may 
only exemplify the crustal field's structure. 
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Fig. 6. Moduli of (s,t)(z) of the fastest growing mode. 
Solid, dash-dotted and dashed lines refer to Bq = 2000, 
Bq = 100 and Bq = 10, respectively. 
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7. Field structure of the fastest growing mode for 
= 2000. Arrows: b X)Z , grey shading: value of b y , dark - 
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An acceptable approximation of the radial profile of a 



dipolar crustal field as given, e.g., in Page et al. (2000) will 
in general have to allow for non-zero values at the bound- 
aries which cause a more complicated perfect conductor 
boundary condition. Moreover, the strong dependence of 
the conductive properties on the radial co-ordinate should 
anyway be taken into account. 

Recently, Hollcrbach fc Riidigcr (2002)) solved the non- 



which support our results: For the largest feasible initial 
field their Legendre spectra are no longer convergent and 
indicate a local minimum at I ~ 60, a hint that small-scale 
field structures are fed in a non-cascade-like, possibly un- 
stable way. Moreover, their findings confirm the impor- 
tance of the background field's profile curvature: For an 
initial (poloidal) field nearly linear in r the excitation of 
small-scale structures and the acceleration of the decay 
were insignificant, in contrast to an initial (toroidal) field 
with a nearly quadratic r-profilc. 

To get an impression of possible consequences for the 
evolution of neutron stars we now simply assume, that the 
real B-profile is sufficiently "curved" (i.e., its second spa- 
tial derivative is big enough) and associate the parameter 
Bq with a typical value of the field. 

Assuming further electric conductivity and chemical 
composition to be constant, gq = 5 x 10 26 s _1 and the rel- 
ative atomic weight A/Z = 25, respectively, we find the 
normalization field at a density p = 10 14 g cm -3 to be 
7 x 10 10 G (see, e.g., [Page ct al. 2000| ). That is, for typical 
(inner) crustal magnetic fields ranging from 7 x 10 12 G to 
1.4 x 10 14 G we find a B between 100 and 2000 and the 
e-folding time of the most rapidly growing unstable mode 
to be 0.0035 and 0.0003 times the ohmic decay time, re- 
spectively. Note, when comparing with the values of Fig. 
that the dimensionless thickness 2 of the slab causes a fac- 
tor 4 in the growth times. These values agree roughly with 
the ratios Ts/r h m derived in Sect. || from observational 
data. An initial perturbation will quickly evolve to a level 
at which the linear analysis is no longer feasible, that is, 
at which it starts to drain a remarkable amount of energy 
out of the background field. Note, that the assumption 
about Bq may well be in agreement with the surface field 
data given in Tabic 0. 

We want to emphasize again that a sufficient curva- 
ture of the background field profile is a necessary condi- 
tion for the occurrence of an unstable behavior. Therefore 
neither the derivation of the well-known helicoidal waves 



linear jilq. ( 22) for the first time in a spherical shell ge- 
ometry, again with the transport coefficients constant in 
space and time, and applying the vacuum and the perfect 
conductor boundary conditions at the surface and at the 
crust bottom, respectively. This is the best model for the 
non-linear magnetic field evolution in the crust of neu- 
tron stars available so far. They made different assump- 
tions about the large-scale initial field (which corresponds 
in some sense to our background field) and report, as all 
the formerly mentioned authors too, about insurmount- 
able numerical problems when the initial magnetization 
parameter, very roughly to be identified with B , exceeds 
« 200. Although their investigations can not reveal the 
existence and the character of the instability presented 
above (too short integration period in comparison with 
the growth time of the unstable modes for the maximum 
initial field considered), they found at least two features 



(whistlers) nor its modification presented in Vainshtcin 



et al. (2000) could reveal it because a homogeneous back- 



ground field was used there. 

With great reserve we may speculate about possible 
observational consequences. For SOINSs, the small-scale 
field modes initially generated or existing in the crust have 
already been decayed and the magnetic field is concen- 
trated almost completely in the large-scale, say, dipolar 
mode. Simultaneously, in the process of cooling the coef- 
ficient m*c/er e , i.e., the normalization field Bn, becomes 



smaller and smaller until the Hall-term in Eq. (22) domi 



nates the linear ohmic term. Of course, as a counteracting 
tendency, the magnetic field decay will continue, increas- 
ingly modified by the Hall-cascade and the gradual onset 
of the Hall-instability. Therefore, we are confronted with 
a competition of one process increasing Bq (via r e ) and 
another diminishing it (via the non-normalized magnetic 
field) . The final decision whether Bq indeed reaches values 
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bation is concentrated in deeper regions of the crust (at 
about 25 % of the thickness). From its pattern one may 
infer shear, both in y- and z- directions, to be the pre- 
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Fig. 8. Heat source density perturbation 2 curl Bo • curlb 
(upper panel) and Lorentz force density perturbation 
curl Bq xb + curl b x Bq (lower panel) of the fastest grow- 
ing mode for B = 2000 in arbitrary units. Upper panel: 
dark and light shading denotes cooling and heating by 
the perturbation, respectively; note the stretched z-scale. 
Lower panel: arrows denote x— and z-force components, 
grey 



vailing type of deformation. Thompson & Duncan (1995) 
discussed a scenario in which the hclicoidal waves, being a 
consequence of the Hall-drift, ignite high energetic bursts 
by cracking the crust. From our results we can hope that 
the unstable modes of the Hall-instability act even more 
efficiently than the at best undamped, but never growing 
helicoidal waves can. 

However, all the questions connected with the actual 
influences of the presented instability on spin-down, cool- 
ing and crust-cracking can be decided only by perform- 
ing non-linear calculations in a realistic spherical-shell ge- 
ometry taking into account a realistic density profile and 
chemical composition of the crust. On the other hand, the 
question at which field strength the action of the Hall- 
instability finally ceases, that is, the field decay is again 
overtaken by ohmic dissipation and, as a consequence, to 
which extent the distribution of observable pulsars in a B— 
P, better, in a B-P-r a diagram is possibly influenced, can 
again only be answered on the basis of realistic non-linear 
magneto-thermal simulations. Whether the scenario of an 
accelerated field decay is in agreement with the observed 
properties of the pulsar population can only be decided 
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0(1000) or even higher has to be left to non-linear coupled 
magneto-thermal calculations. Once the Hall-term dom- 
inates the instability may raise small-scale perturbations 
down to scale lengths > Z cr it on expense of the dipolar 
mode, resulting in a spin-down behaviour of SOINSs as 
discussed in Sect. ||. 

Another possible observational consequence follows 
from the dissipative effects of the rapidly excited small- 
scale field modes, that is, enhanced Joule heating. As can 
be seen in the upper panel of Fig. ||, the growing perturba- 
tions cause heat sources located close to the surface of the 
crust which, in turn, may produce a patchy neutron star 
surface. Note, that the plot shows only the part linear in 
b of the heat source distribution, since the quadratic term 
(curlb) 2 has to be dropped in the linear analysis. These 
heat sources are concentrated in a very thin layer com- 
prising only less than 5 % of the crustal thickness. Clearly, 
conclusions about the observability of that hot spots are 
impossible to be drawn from our simple model, since to 
infer the temperature distribution from the heat sources 
demands the knowledge of the absolute perturbation am- 
plitudes. Anyway, we expect the neutron star's surface to 
be warmer during the episode of the Hall-instability than 
standard cooling calculations predict. 

Third, strong small-scale field components cause 
small-scale Lorentz forces. The lower panel of Fig. || shows 
that the maximum of the Lorentz force density pertur- 



by help of a population synthesis as presented by Rcgim 



bau & dc Frcitas Pacheco (2001) applying a decay pattern 



with an initially slow decay followed by an episode of rapid 
decay and leading finally to a nearly constant field. 
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